Stellar Dynamics and Black Holes 

DAVID MERRITT^ 

^Rochester Institute of Technology, Rochester, NY, USA 



1. Introduction 

Chandrasekhar's most enduring contribution to stellar dynamics is probably dynamical 
friction. The history of that discovery is fascinating and well known |17]. The first 
(1942) edition of Principles of Stellar Dynamics contained a chapter entitled "The Time 
of Relaxation of a Stellar System," in which Chandrasekhar presented his "completely 
rigorous" calculation of the rate at which a star is accelerated by encounters with other 
stars. The result was expressed in terms of ^ AE'^, the accumulated squared kinetic 
energy change after some time At. From this, Chandrasekhar defined the relaxation time 
as 

J2AE'^ _ At 

and showed that Te was very long, of order 10^^ yr, for a system like the Milky Way. He 
then defined a second relaxation time , Td, in terms of the angular deflection of an orbit 
and showed that it was equivalent to Te- 

It was apparently only after publication of his monograph that Chandrasekhar recog- 
nized another important consequence of gravitational encounters. In the absence of a de- 
celerating force, repeated encounters would increase the kinetic energy of a star without 
limit. In a fluid medium. Stokes' law states that there is a frictional force, proportional 
to the velocity, that competes with the acceleration. In a paper published in 1943 fl, 
Chandrasekhar recognized the necessity for such a force in the stellar dynamical case, 
and showed that the assumption of a Maxwellian velocity distribution implied a simple 
relation between the frictional force and his previously-derived coefficient of diffusion. 
He then derived the coefficient of dynamical friction by summing the velocity changes 
experienced by a test body, in the direction of its motion, as field stars were deflected 
around it. 

Chandrasekhar was clearly proud of the fact that Brownian motion in stellar systems 
could be treated in such an exact way, compared with the "intuitive character of the as- 
sumptions" in the case of fluids He wrote: 

The discussion of the physical foundations of the theory of Brownian 
motion [in fluids] in the preceding sections has disclosed certain inherent 
limitations in the theory. The limitations are nowhere more serious than in 
the circumstance that the coefficients q and 77 [the coefficients of diffusion 
and friction] are not derived from a microscopic analysis of the individual 
encounters. It is therefore of interest that stellar dynamics provides a case of 
Brownian motion in which all phases of the problem can be explcitly ana- 
lyzed. 

In this lecture, I briefly review Chandrasekhar's derivation of the dynamical friction 
cofficient, then discuss its application in a few cases that are relevant to the nuclei of 
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galaxies containing massive black holes. I conclude by discussing recent work that ex- 
tends Chandrasekhar's ideas to star clusters in which the effects of relativity can not be 
ignored. 



2. Dynamical Friction 

In deriving the velocity diffusion coefficients in a stellar system, Chandrasekhar assumed 
that the instantaneous forces acting on a test star could be separated into two components: 
a component due to the smoothed-out distribution of matter, and a component arising 
from chance stellar encounters. He assumed that the encounters were independent, and 
computed their net effect by integrating the momentum changes of single encounters over 
the variables defining the relative orbit of the binary (test star, field star) system. The 
smoothed-out component was assumed to be infinite and homogeneous, implying a con- 
stant velocity for the test mass in the absence of encounters. Aside from enforcing istropy 
in velocity space, no restrictions were placed on /(tt), the distribution of field star veloc- 
ities at infinity. 

Integrating over all impact parameters p, the change per unit interval of time of the test 
star's velocity, in a direction parallel to the initial relative motion, is then given by 



2TTG^m f (to / + to) n _ 



(Az;||) = ^ In {l+piUPo) ■ (D 

Here to and irif sae the masses of the test and field star respectively, n is the field star 
density, V = v — uis the relative velocity of test and field star before the encounter, and 

G (TOf + to) 

P0 = 

is the impact parameter corresponding to an angular deflection of 7r/2. The final step is 
the integration over field star velocities. Since equation ([T]i gives the velocity change in 
the direction of the initial relative motion, we must multiply it by 

V ■ V V — Ux 



Vv V 

to convert it into a velocity change in the direction of the test star's motion, assumed here 
to be along the x axis. The dynamical friction coefficient is then 



■ V — Ux 



(Az;||) = / /(„)(Am||)— 



{mf+m)p J /(«)!i-^ln 



1 



„2 yi 

ymax ^ 



{mj + rnf 



du 



where p = mfU. If the field star distribution is isotropic in velocity space, there are 
various ways to simplify the integral over velocities. Chandrasekhar chose one that was 
especially propitious: representing the velocity-space volume element in terms of u and 
V yields 

167r^G^(TO;+TO)p r . 2.,, w 
(Au||) = — / j(u)u H (u, u,p,„ax)au, 

2^ I 1 max ^ 







G-^{'mf + to)2 



(2) 
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This choice separates the integrand into a product of / with a weighting function H (which 
Chandrasekhar wrote as J/8u). 

Chandrasekhar then pointed out an approximation that greatly simplies "Under 
most conditions of practical interest," he wrote, the argument of the logarithm in equa- 
tion (|2]i "is generally very large compared with unity." Invoking this assumption, Chan- 
drasekhar derived various approximate forms for H. The last, and simplest, of these was 



n 



In A if u > M, 
if V < u 



which corresponds to treating the logarithmic term in equation (|2| as a constant. In this 
approximation, the dynamical friction coefficient becomes 

IGtt^G^ {mf + m)plnA r ^ 
(Awii) = / f[u)u du. (3) 







Equation Q contains the beautiful result that the dynamical friction force is proportional 
to the number of field stars moving more slowly (at infinity) than the test star. 
For the "Coulomb logarithm", Chandrasekhar suggested the form 



InA = In 



G{mf + to) 



As for the parameter Pmax, Chandrasekhar and von Neumann H argued that the effect 
of gravitational interactions should be cut off at roughly the interparticle distance. Later 
authors [e.g. 7] have generally advocated larger values for Pmax^ the physical size of the 
system, or the density scale length. 

Following the derivation of equation (O, Chandrasekhar wrote: "As we shall see presently, 
it is precisely on this account [i.e. on account of the fact that only stars with u < v con- 
tribute to the force] that dynamical friction appears in our present analysis." 

Maybe I missed it, but I could not find where in his paper Chandrasekhar justifies this 
statement. In any case, it is fun to play devil's advocate, and ask how dynamical friction 
acts in cases where this is not true. 

A student of Chandrasekhar's, Marvin Lee White, considered one such case ll20ll . If 
Pmax is equated with the half-mass radius of a stellar cluster, then using the Virial Theo- 
rem, it is easy to show that InA w In(A^) with N the number of stars in the cluster For 
sufficiently small N, the assumption of large In A breaks down. White showed that in 
a cluster containing just a few hundred stars, a significant fraction of the frictional force 
acting on a test mass could come from stars with u > v. 

Two other cases come to mind, both in the context of massive black holes (MBHs). 



I A MBH near the center of a galaxy experiences a kind of random walk due to 
perturbations from stars. Assuming equipartition of kinetic energy, its velocity 
will be small, v « (TOi/M)^/^\/u2 jf y small, the logarithmic term in 
equation (|2]l will be close to zero for all field stars with u < v. Either the frictional 
force acting on the MBH is much less than implied by equation (|3]l, or most of the 
force must come from field stars with u > v. 



The integral for "H can be expressed exactly in terms of simple functions, although the expression is quite 
long. This is much easier to demonstrate nowadays than it was in the 1940s, with tools like Mathematica. 
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2 There are physically reasonable models for the distribution of stars around a MBH 
that imply few or no stars with velocities less than some value; for instance, the 
local circular velocity. In such cases, Chandrasekhar's approximate formula would 
predict no frictional force on a body moving in a circular orbit; all the force would 
have to come from the fast-moving stars. 

Consider first the Brownian motion example. Chandrasekhar (1943) showed that the 
rms velocity of the test star is given by 



3 

V — — 

rms 2 



((A« 



{^y\\)ly 

Expanding the diffusion coefficients about v — Q, 

(Au||) = -Av + Bv'^ 
((Aw||)2) ^ C + Dv\.. 

so that w^jjjg — 3C/2A. The coefficients A and C are obtained from the low-velocity limit 
of equation (|2]l, and the equivalent expression for ((Aw|| )'^) The results are 

Ifi f 11^ 







C = -n'G'm^p / f{u) log ( 1 + ) u du (4b) 




where AI is the MBH mass. It is evident from the first of these expressions that field 
stars of every velocity contribute to the dynamical friction force. If the field star velocity 
distribution is a Maxwellian, the result of the integrations is 

A = 4V6^G2Mp In A', C = ^V^G'^m^p ^ In A' 

where 

In A' = — / dz e , ^ , „„ ^ , 

2 io \ 9G2Af2 / V 9G2M2 

The rms velocity of the MBH has the expected value, (rrii^u'^ /Aiy^^, regardless of the 
value of In A' . But returning to the question of which stars generate the frictional force: 
the integrand of equation ( l4al i peaks at u « (GAf/pmax)^^^- Furthermore, Pniax in the 
case of a MBH near the center of a galaxy is expected to be comparable to the MBH 
influence radius GM/u^ 1 13, 15]. It follows that most of the frictional force responsible 
for Brownian motion of a MBH comes from field stars with speeds similar to their rms 
value, and not exclusively, or even predominantly, from the slow-moving stars (Figure 1). 

The second case mentioned above was a test mass orbiting around a MBH at the center 
of a galaxy. Suppose that the stars, which are responsible for the frictional force, are 
distributed around the MBH with density p{r) oc r~'' . Assuming an isotropic velocity 
distribution, their phase space density is given by 

fo,{~Efr'-^l^o,{2vl-uy" (5) 

where Vc = (GMlr)^!"^ is the local circular speed. For 7 < 3/2, / diverges at u = 
2^/2^^^ and in the limit 7 = 1/2, / is a delta-function at u = 2^/2^^ other words, in 
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Figure 1. Relative contribution of different field-star velocities to the dynamical fric- 
tion integral, equation Bat , in the limit of low and high velocities of the test mass, 
assuming a Maxwellian distribution of field-star velocities with one-dimensional ve- 
locity dispersion cr/ = (u^/S)^^^. pmax = GAI/aj has been assumed, appropriate 
for a massive test body moving at the center of a galaxy. Roughly the same set of field 
stars dominates the frictional force, whether the test mass is moving very slowly or very 
fast. 



a p cx r^^/^ density cusp around a MBH, none of the stars at r is moving more slowly 
than the local circular speed. The dynamical friction formula in its approximate form (|3) 
would predict zero frictional force in this case. 

It so happens that the distribution of stars at the centers of galaxies containing MBHs 
is often as flat as p ~ r~^/^ |@]; indeed this appears to be the case at the center of the 
Milky Way [jj. One explanation is that these low-density cores were "carved out" by 
binary MBHs ifldilFI But whatever their origin, the argument just given suggests that 
dynamical friction in such cores must be due almost entirely to stars that move faster than 
the inspiralling body. Of course, the frictional force on a test mass at a distance r from the 
MBH would come partly from stars at greater r, and some of these stars will have smaller 
u. Chandrasekhar's theory does not make clear predictions in the case of inhomogeneous 
systems. It would be interesting to address this question using an A^-body code. 

Why do the two populations - stars with u < v and stars with u > v - contribute in 
such different ways to the frictional force? One way to visualize this is to compute the 
dynamical-friction "wake," the overdensity of stars behind the test body that is responsible 
for the decelerating force. If the field star distribution is infinite and homogeneous, and if 
the test body's velocity is rectilinear and unaccelerated (all assumptions that were made 
by Chandrasekhar), the wake is time-independent in a frame following the test mass, 
and its density can be computed using Jeans's theorem given the distribution of field star 
velocities at infinity lfl6ll . 

Figure 2 shows the result of such a calculation, assuming that the velocity distribution 

^If so, then one might expect / to be somewhat ajiisotropic. 
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Figure 2. Dynamical friction wake around a massive object moving tlirougii an infinite 
homogeneous stellar system. The panels at the top show contours of the density, in a 
frame that follows the massive body (located at the origin); motion is to the right at con- 
stant speed V. The left panel shows the total density, the middle panel shows the density 
contributed only by stars with velocity at infinity less than that of the massive body, and 
the right panel shows the contribution from the complement of stars that move faster 
than V. The lower panels show the density along a line through the moving body in the 
direction of its motion. In the usual approximation, first derived by Chandrasekhar, one 
computes the dynamical friction force ignoring the fast-moving stars. This figure helps 
to explain why that is a good approximation, even if the fast stars dominate the total 
density. 



at infinity is given by equation (|5]l with 7 = 5/4 and that the test body's velocity is Vc- 
For this choice of /, the "fast" stars dominate both the total density at infinity, as well as 
the density in the wake. However the shapes of the two, partial density wakes are very 
different. The wake created by the fast stars is elongated counter to the direction of the 
test body's motion, reducing the net force that it exerts on the test mass. In addition, the 
change in density between the upstream and downstream sides of the test mass is less for 
the fast stars than for the slow stars. These two differences are responsible for the small 
contribution of the fast stars to the total frictional force, in spite of their higher density at 
infinity and in the wake. 



3. Relativistic Dynamical Friction 

In 1969, Edward E. Lee, a PhD student of Chandrasekhar's, re-derived the expression for 
the dynamical friction force, using the first post-Newtonian (IPN) approximation for the 
relative orbit of test and field stars ||I2|| . Lee motivated his calculation in the following 
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way: "These effects of general relativity, in the lowest order, may be relevant for consid- 
eration of stellar encounters in dense systems such as are now contemplated in a number 
of contexts." Lee's final expression for the relativistic dynamical friction coefficient was 
fairly complex, and he did not go so far as to evaluate it in the case of a particular velocity 
distribution. But given f{u), Lee noted that the relativistic formula predicted a higher 
frictional force than the non-relativistic formula, due to the greater relative deflection in 
the IPN approximation, and also because the relativistic transformation from the inertial 
frame to the frame of the test mass introduces an additional term into the force. 



4. Relativistic Stellar Dynamics 

Massive black holes were not yet quite in vogue at the time that Lee wrote his thesis. It 
is now commonly assumed that all massive galaxies have them. Given Chandrasekhar's 
fundamental contributions to the theory of gravitational encounters, and to the theory of 
black holes [5], gravitational encounters between stars orbiting near a MBH is a topic that 
he might naturally have been interested in. 

A black hole is a compact object, and from the point of view of a star that orbits far 
outside the event horizon, its gravitational field should be nearly indistinguishable from 
that of a Newtonian point mass. But the Newtonian approximation breaks down for matter 
that passes within a few gravitational radii: 

Tg EE — - « 5 X 10-*pc 



since the orbital velocity at such distances approaches the speed of light. 

At first sight, one is struck by the enormous difference between rg and r^afi, the radius 
of gravitational of influence of a MBH: 

GM f M \ f 

nnfl = — 5- « Ipc ' 



\W^Mq) V 100km 

where o-^ is the stellar velocity dispersion at r > rjnfi. These two equations suggest 
that the vast majority of stars within the influence sphere are too far from the MBH for 
relativity to be important. 

This conclusion turns out to be misleading, for a couple of reasons. First: the effects 
of relativity depend less on the size of an orbit than on its distance of closest approach 
to the MBH. The lowest-order corrections to the Newtonian equations of motion have 
amplitudes that are of order V^^ where V is the penetration parameter. 

Here, a and e are the semi-major axis and eccentricity of the orbit, and rp = (1 — e)a 
is the radius of periapsis. It turns out that the feeding of stars and compact objects to 
MBHs occurs predominantly from very eccentric orbits. For instance, capture of stellar- 
mass BHs by MBHs - so-called "extreme-mass-ratio inspiral" (EMRJ) - is believed to 
take place from orbits with semi-major axes a < 0.1 pc and eccentricities in the range 
0.999 — 0.99999 [ H]. For such orbits, V can be of order unity even though the orbit 
extends outward to thousands or tens of thousands of gravitational radii. 

A second reason why the effects of relativity can not be ignored has to do with the way 
in which stars get placed onto orbits of such high eccentricity. The dominant mechanism 
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is believed to be torques, i.e. non-radial forces, that arise from the slightly aspherical 
distribution of matter near a MBH fil8il . These torques remain effective as long as orbits 
near the MBH - both the orbit of the star being torqued, and the orbits of the torquing stars 
- maintain their orientations; any mechanism that causes orbits to precess (for instance) 
tends to randomize the torques. Relativistic precession of the periapsis - or, as it is know 
in the Solar system, precession of the perihelion - is such a mechanism. If the time scale 
for relativistic precession is shorter than the time time required for the torques to do their 
work, feeding of objects to the MBH will be greatly inhibited. This relativistic quenching 
effect turns out to be of major importance in the EMRI problem. 

Consider a star on a bound orbit near a MBH, with semi-major axis a and angular 
momentum L. The orbit-averaged torque | x r | that another star of mass m exerts on it is 
~ Gm/a. The residual torque due to N stars, randomly oriented about the MBH at about 
the same radius, is \T\ w ^/N Gm/a. Over some span of time, this orbit-averaged torque 
is nearly constant, and the angular momentum of a test star's orbit responds by changing 
linearly with time, AL/ At w T w ^/NGm/a. This coherent evolution continues for a 
time ^ tcoh, where tcoh is the time scale associated with the most rapid mechanism that 
randomizes the torques, i.e. the orientations of the orbits resposible for the torques. 

Sufficiently far from a MBH, coherence breaking is due mostly to the same distributed 
mass that generates the torques. Modelling that mass as spherical, the associated preces- 
sion time is 



where P is the orbital period, AI is the MBH mass, and A/* is the distributed mass within 
r — a. (There is also a dependence on e, in the sense that eccentric orbits precess more 
slowly.) This time is believed to be of order 10"' yr in the case of the so-called S stars, 
bright young stars at the center of the Milky Way whose orbits can be tracked astrometri- 
cally; the S stars have 0.01 pc < a < 0.1 pc and lO^r <P< lO^yr fl. 

Closer to the MBH, relativistic effects become important. The most important preces- 
sional time scales associated with relativistic corrections to the equations of motion are 



The first of these refers to precession of the periapse; the subscript S indicates that this 
effect is due to the Schwarzschild (zero spin) part of the metric. This "Schwarzschild 
precession" is similar to, but in the opposite sense of, precession due to the distributed 
mass. The other two time scales are associated with the spin (J) and quadrupole moment 
(Q) of the MBH. tj is the time for precession of the line of nodes due to frame-dragging, 
and tq is the nodal precession time resulting from the nonsphericity of space-time around 
a spinning hole. (The quadrupole precession rate also depends on the inclination of the 
orbit with respect to the MBH spin axis.) These latter two times are functions of x, the 
dimensionless spin of the MBH: writing the hole's spin angular momentum as J, 




(6) 




tj 



GAP 
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and X = 1 for a maximally-spinning hole. 

The time scales for relativistic precession defined above decrease as (1 — e), or faster, 
as e — > 1. If torques from the ^/N asymmetries in the stellar distribution drive the eccen- 
tricity of a test star's orbit to a sufficiently large value, relativistic effects will dominate 
its precession. At some critical precession rate (i.e. eccentricity), the sign of the torque as 
experienced by the orbit will fluctuate with such a high frequency that its net effect over 
one precessional period will be negligible: in other words, relativity will "quench" the 
effects of the torques. 

This critical eccentricity can be estimated as follows. As shown above, the residual 
torque produced by N stars is 

a \/N{a) a 

where N{a) = Mi,{a)/m is the number of stars, of mass m, within radius r = a. Writing 

L = [GMa{l — e^)] for the test star's angular momentum, the time scale over which 
L is changed by this torque is 



1 dL 



M 



GM 



The condition that this time be shorter than the relativistic precession time ts is 

r„ M 



a M^(aj 

where 

i = L/U = (1 - e2)i/2 

is the ratio of the test star's angular momentum to its circular value. In terms of eccentric- 
ity, 

2 /mV 1 



2 X 10"^ 



a \ ( M \ ( m \ ( N 



o.oipcy \iq^MqJ \ioMqJ \io 



^2 



This "Schwarzschild barrier" (so called since it arises from the Schwarzschild part of the 

metric) sets an effective upper limit to the eccentricty of the test star. 

What happens when a star (or stellar remnant) "strikes" the barrier? Figure |3] shows 

I — I 

one example, extracted from an A^-body simulation |14]. The orbit's eccentricity first 
oscillates, about a lower bound given by equation (O. The oscillations have a period 
equal to ig; they reflect the periodically changing effects of the torques on the precessing 
star After some elapsed time - roughly 6 precessional periods in this case - the star is 
"reflected" from the barrier by the torques, back to smaller values of the eccentricity. The 
elapsed time until reflection is just the coherence time tcoh defined above. The orbits of 
all the other stars are also precessing: not as rapidly as the test star (which has the highest 
eccentricity), but at a lower rate - in the example shown in Figure |3] the distributed mass 
determines the precession rate of most stars, via equation (|6]l. After a time of tcoh ~ ^m, 
the net torque due to the other stars has changed its direction and the angular momentum 
of the test star responds by decreasing, taking the star away from the barrier. 
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semi-major axis argument of periapsis log (1 - e) 

10 ^^-^ — I 1^ — ^ — -I 1— .— ^ — I I I . 




1.8 2 2.2 1.8 2 2.2 1.8 2 2.2 



time 

Figure 3. An eccentric orbit near a massive black hole, extracted from an A'^-body in- 
tegration [14]. Plotted are the semi-major axis, argument of periapse and eccentricity 
versus time, at low (top) and high (bottom) time resolutions. In the plots of eccentricity 
vs. time, the lower (red) curves show the predicted location of the Schwarzschild bar- 
rier. Changes in the predicted barrier location reflect changes in the semi-major axis, 
via eq. jSjl. 
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The iV-body experiments from which Figure |3] was derived were designed to generate 
EMRIs, i.e. inspiral of stellar remnants into the MBH. For this to happen, the remnants 
must sometimes penetrate the Schwarzschild barrier 

Penetration does occur, but because the coherent torquing - or "resonant relaxation" 
- is quenched by the effects of relativity, all that is left is the kind of random gravita- 
tional scattering treated by Chandrasekhar, or "non-resonant relaxation". It turns out that 
the rate at which objects get scattered past the angular momentum barrier is accurately 
predicted by Chandrasekhar' s formulae. The only additional consideration arises from 
the fact that stars remain near the barrier only for a limited time. At w tcoh, before the 
changing background torques pull them away; to be effective, the encounters described 
by Chandrasekhar' s formulae must decrease the test body's angular momentum, by an 
amount greater than the amplitude of the oscillations shown in Figure |3] in a time less 
than tcoh- This argument correctly reproduces the EMRI capture rate observed in the 
simulations, and also turns out to imply a minimum semi-major axis, of order 10"'^ pc, 
below which stars can not be driven past the barrier, either by resonant or non-resonant 
relaxation 

The collisional dynamics of relativistic stellar systems is a fascinatingly rich topic. It 
is a shame that Chandrasekhar is not here to sort it out for us. 
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